Shear-induced fractures and three-dimensional motions in an organogel 
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The flow behavior of a viscoelastic organogel is investigated using ultrasonic velocimetry combined 
with rheometry. Our gel presents a decreasing flow curve, i.e., the measured stress decreases as a 
function of the applied shear rate. Strikingly, we note that the local flow curve calculated from 
the velocity profiles also exhibits a decreasing part. We attribute this regime to the presence of a 
fracturing process and three-dimensional motions in the bulk of the sample. 
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I. INTRODUCTION 

Soft glassy materials are disordered materials with very 
long relaxation times. They encompass many differ- 
ent systems such as foams, emulsions, granular systems, 
pastes, and gels. The response of such systems to an ex- 
ternal shear stress is characterized by two regimes: for 
stresses below the yield stress ctq they lose their ability 
to flow, remain jammed, and respond elastically, whereas 
for stresses above ctq they flow as liquids [l| . Despite this 
macroscopic shared behavior, local velocimetry experi- 
ments have revealed that the flow of jammed materials 
is not universal Nonadhesivc emulsions micro- 
gels, and pastes were shown to flow homogeneously in 
the vicinity of the yielding transition [jl. In contrast, 
Laponitc gels 0| or adhesive emulsions 01 display shear 
localization as they go through the yielding transition Q . 

In some particular cases, the applied shear stress de- 
creases as a function of the applied shear rate and the 
flow curve presents a minimum 0, Q. In a pioneering 
work, Pignon et al. Q have studied thixotropic colloidal 
suspensions displaying such a behavior. They have mea- 
sured the local strain by introducing a fine filament of 
colored suspension in the sample and following its defor- 
mation. Four flow regimes corresponding to four differ- 
ent states of the strain field were defined by combining 
visualization techniques and rheometric measurements in 
cone-and-plate geometry. At small deformations, the gels 
have an elastic behavior corresponding to a homogeneous 
strain field. Above a critical strain, the gel is deformed 
as a function of the applied shear rate. At very low shear 



rates, the thickness of the sheared layer is of the order 
of the size of the clay particles and the shear stress de- 
creases as a function of the apparent shear rate. In the 
intermediate shear rate range, the size of the fluid layer 
increases and the shear stress remains constant as the ap- 
parent shear rate increases. At high shear rates, the flow 
becomes homogeneous. This points out that the global 
flow curve does not reflect the local rheological behavior 
and raises the question of the origin of decreasing flow 
curves. 
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In this work, we revisit this problem using an organogel 
and ultrasonic velocimetry in Couette geometry ^ com- 
bined with standard rheometry. We show that three- 
dimensional motions and fractures occur in the decreas- 
ing part of the flow curve. This limits the possibility 
to etablish a valid flow curve even by using local tools 
such as local velocimetry. Indeed, it is not possible to in- 
duce pure shear flow (even locally) in this system. This 
paper is organized as follows. Section II introduces the 
organogel under study, a mixture of an organic solvent 
and a low molecular weight gelator. We describe the syn- 
thesis of this system and study its structure as well as its 
basic linear rheological properties. In Section III, the 
flow properties of this thixotropic gel are addressed. We 
first use nonlinear rhcology alone to access the flow curve 
of our system. Then the ultrasonic setup for measuring 
velocity profiles under shear is briefly recalled. Finally, 
evidences for fractures, three-dimensional motions, and 
a multivalued, decreasing local flow curve are provided 
and discussed. 
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FIG. 1: The N-5-hydroxypentyl-undecanamide molecule. 



II. THE ORGANOGEL UNDER STUDY: 
SYNTHESIS, STRUCTURE, AND LINEAR 
RHEOLOGY 

The present work is devoted to the study of an 
"organogel." Organogels are thermoreversible, viscoe- 
leastic materials consisting of low molecular weight gela- 
tors (LMWG) self-assembled in an organic solvent into 
complex three-dimensional structures. At high tempera- 
ture, LMWG are quite soluble in organic solvents and 
their solutions are liquid [HI. At low tempera- 
ture, multiple nonconvalent interactions such as hydro- 
gen bonding, donor-acceptor, and hydrophobic interac- 
tions between the organogel building blocks give birth to 
a three-dimensional network of entangled and connected 
fibers. This process is reversible, i.e the gel can be melt 
again at high temperature and reformed under cooling. 

In recent years, these organogels have been the sub- 
ject of increasing interest mainly due to the fundamen- 
tal questions raised by the gel structures and assembly 
mechanisms HI, E El, [H 13. In this paper, we fo- 
cus on their rheological properties. A mixture of N- 
5-hydroxypentyl-undecanamide (N5HU) in toluene has 
been studied. We first describe briefly the synthesis of 
N5HU and the gel formulation. Then the gel structure is 
characterized together with its linear rheological proper- 
ties. 



A. Synthesis of N-5-hydroxypentyl-undecanamide 

The gelator molecule used in this work, namely N- 
5-hydroxypentyl-undecanamide (N5IIU) is shown in fig- 
ure [TJ The synthesis of N5IIU was previously described 
in Ref. [13] ■ Lauryl chloride (1 mol equiv) is added drop- 
wise at 0°C to a solution of 5-aminopentanol (Sigma) (1 
mol equiv) in 40 mL of 0.2 N NaOH (Aldrich). The re- 
action is stirred for one night at room temperature. The 
precipitate is filtered and dissolved in dichloromethane 
(Aldrich) . The organic layer is washed several times with 
water and dried over Na2S04 (Sigma). After evapora- 
tion, the white solid is triturated with pentane, filtered, 
and dried under a vacuum. Toluene (Aldrich) is used as 
received. Purity of the sample was checked by IH NMR 
and IR spectroscopy. 

This molecule is derived from previous studies on N- 
acyl-1, w-amino acid based gelators [l^ and on N-3- 
hydroxypcntyl-undecanamide [l9|. It forms a gel in var- 
ious fluids such as toluene and some mixtures of silicone 



60 
50 
40 

y 30 

20 
10 




Liquid 



Gel 



Suspension 



2 3 4 

C wt% 



FIG. 2: Temperature-composition phase diagram for the 
NSHU-toluene system. 



oil in toluene and of pcntanol in dodecane. In the follow- 
ing, we present a study of gels of N5IIU in toluene. 



B. Gel formulation and phase diagram 

To obtain a gel, heterogenous mixtures of various con- 
centrations of gelator N5HU (from 1 to 4.7 wt %) in 
toluene were heated (at about 70° C) until clear liquids 
were obtained. The homogenous warmed solutions were 
then cooled rapidly at a given temperature thanks to a 
temperature-controlled bath. 

Figure [5] shows the phase diagram obtained by study- 
ing the transition temperature between the liquid and 
gel phase for the different samples. The procedure used 
to determine the transition temperature is rather crude. 
The solution is prepared in a tube which is held upside 
down once the final temperature is reached. If the sam- 
ple flows, it is called a liquid, whereas if it sustains its 
own weight, it is called a gel. For concentrations above 
1.2 wt %, gels are obtained at low temperature and be- 
come liquid at high temperature. The gel can be melt 
again at high temperature and reformed under cooling, 
so that the gclification process is reversible. The melt- 
ing temperature and the gclification temperature differ 
by less than 1°C. For concentrations lower than 1.2 wt % 
of N5HU; a suspension is observed at low temperature. 

In order to characterize the structure of the gel and 
to understand why a gel is observed at low temperature, 
microscopic observations and linear rheology experiments 
have been performed and are presented in the next sec- 
tions. 



FIG. 3: Optical microscopy picture of a 1.2 wt % N5HU- 
toluene gel formed at a quench temperature of 5°C. The pic- 
ture was taken just after the quench. The scale bar represents 
100 fim. 



FIG. 4: Optical microscopy picture of a 1.2 wt % N5HU- 
toluene gel at 5°C one week after the quench. The scale bar 
represents 100 /im. 



C. Microscopic observations 

A solution of N5HU at 1.2 wt % in toluene is heated 
at 65°C and cooled down to 5°C at a cooling rate of 
20°C / min. Microscopic pictures of the resulting gel taken 
immediately after the quench and one week later are dis- 
played on Figs. [3] and [4] respectively. 

This qualitative study leads to the following conclu- 
sions. The gel is formed of large aggregates composed of 
interconnected branches or fibers that have grown from 
the same center. The size of the fibers ranges from 20 to 
200 ^m. These aggregates form a network that is respon- 
sible for the elastic behavior of the gel. Figure |4] shows 
that the initial gel state is metastable and that the final 
state (after 1 week) is a suspension of independent and 
unconnected rods which has lost its clastic properties. 
With time some parts of the primary gel fibers melt and 
reform at other places with a different morphology. The 
mean diameter of the fibers increases, the new stucturc is 
much less branched, and Oswald rippening seems to oc- 
cur and to smooth out the stucture. This ageing behavior 
is similar to the one observed by Lescanne and coworkers 
on N-3-hydroxypentyl-undecanamide-toluene gels [3 ■ 

Such an evolution is very slow, which allows us to con- 
sider that the structure does not vary significantly in 
the first five hours. All our experiments will deal with 
fresh samples and will be started a few minutes after the 
quench. Note that in the following, the rheological exper- 
iments will be performed on gels with higher N5HU con- 
centrations. We shall assume that the structure involved 
in these systems is similar to the one presented above for 
lower concentrations. Indeed, due to the turbidity of the 
sample, we have been unable to obtain optical pictures 
for more concentrated systems. 



D. Linear rheology 

To better characterize the nature of the intercon- 
nection between the aggregates, rheological experiments 
have been performed using a stress-controlled rhcometcr 
(TA Instruments AR lOOON) equipped with a cone-and- 
plate geometry made of aluminum (40 mm diameter, 4° 
cone angle, and 56 ^m gap truncation). To avoid sol- 
vent evaporation, a solvent trap was placed above the 
cone. A striated cone-and-plate geometry was used to 
minimize wall slip. A heated solution (at around 65°C) 
of 4.7 wt % of N5IIU in toluene is loaded in the shear- 
ing gap and cooled down to 5°C using the Peltier plate 
of the rheometer. A gel is obtained and the rheological 
experiment is started 10 minutes after loading the sam- 
ple. In order to define the so-called linear regime, the 
sample is first submitted to a periodic stress oscillating 
at a frequency of 1 Hz and with an amplitude a that 
increases exponentially with time. During this experi- 
ment, the storage modulus G' and the loss modulus G" 
are recorded as a function of the stress amplitude and 
plotted in Fig. [5l 

At low stress values, i.e., in the linear regime, both 
dynamic moduli arc independent of the stress amplitude 
a and reveal the properties of the unpertubed gel: G" is 
about 10 times larger than G" which reflects the domi- 
nant character of the material. Up to large experimental 
scatter on G" below 1 Pa, both moduli remain roughly 
constant below a critical shear stress of about 3 Pa, which 
can be taken as the limit of the linear domain. Between 
CT = 3 Pa and 60 Pa, G' slightly decreases as a function 
of the applied shear stress. Above u = 60 Pa, G' drops 
dramatically while G" decreases sharply. This behavior 
is due to the partial breakup of the gel. Thus our N5HU 
gel behaves as a weak gel [201 ■ 

In a second set of experiments, G' and G" are mea- 
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FIG. 5; Viscoelastic moduli G' (•) and G" (o) as a function 
of the amplitude a of the applied shear stress for a 4.7 wt % 
N5HU-toluene gel. The frequency of the stress oscillations is 
1 Hz. 



FIG. 6: Frequency sweep experiments for a 4.7 wt % N5HU- 
toluene gel: G' (• and T) and G" (o and o) as a function of 
frequency. The amplitude of the applied shear stress is 2 Pa, 
which falls into the linear regime. 



sured as a function of the angular frequency (which is 
varied from 0.1 to 10 Hz) for a fresh sample of the N5HU 
toluene gel and a stress amplitude of 2 Pa. This stress 
falls into the linear regime so that the frequency sweep 
is nondestructive and the gel structure remains intact. 
Two different measurements performed on the same sam- 
ple within 20 minutes are reported in Fig. [H] showing a 
good reproducibility of the experiment. A slight increase 
of the clastic modulus G" as a function of frequency is 
noticed (from 1.5 lO'' Pa at 0.1 Hz to 2 10'^ Pa at 10 Hz), 
which is consistent with a viscoelastic behavior. On the 
other hand, the loss modulus is observed to increase as 
the frequency decreases. Since G" should reach zero at 
zero frequency, a viscoelastic process occurs with a very 
long characteristic time, corresponding to a very low fre- 
quency value (< 0.1 Hz). This long relaxation time may 
be attributed to local rearrangements of the aggregates. 

To conclude this study, it is important to note that 
reproducibility is excellent on the same sample in the 
same day (as shown by Fig. [6]) but not as perfect when 
we compare the properties of two different gel samples 
loaded and formed in the rhcomcter. The values of G' 
and G" may vary by 40% from one experiment to an- 
other. This behavior is common to many organogels. 
We believe that these discrepancies are related to a dis- 
persion in the aggregate number density from one gel 
to another. Indeed, this parameter is very sensitive to 
cooling velocities, impurities, and filling conditions of the 
geometry which arc difficult to control very precisely. Af- 
ter studying the main linear rheological properties of our 
organogel, we now focus on its nonlinear rhcology. 



III. NONLINEAR FLOW BEHAVIOR: 
EVIDENCE FOR FRACTURES IN THE BULK 
AND THREE-DIMENSIONAL MOTIONS 

In this section, we study the nonlinear rheological re- 
sponse of our NH5U-toluene gels. We first address the 
question of the sensitivity of these gels to shear. Indeed, 
most organogels are very sensitive to shear stress which 
tends to expel the solvent from the three-dimensional net- 
work and to break the gel. When the mechanical excita- 
tion is stopped, the sample behaves as a suspension and 
has lost its clastic properties. Usually, the only way to 
recover a gel is to heat it and cool it down once again. 
Here we first present data showing that this is not the 
case for our NH5U-toluene gels and that recovery occurs 
spontaneously once shear is stopped. After discussing 
global nonlinear rheological properties, we present local 
velocity measurements using ultrasonic velocimetry and 
extract the local rheological flow curve. 



A. Sensitivity of the gel to shear stress 

To study the flow properties of the gel and its recov- 
ery, we go back to the experiment presented in Fig. [5] 
and performed on a 4.7 wt % NH5U-toluene gel in the 
striated cone-and-plate geometry. Figure [7] presents the 
storage modulus G' and the loss modulus G" measured as 
a function of time. During the first 480 s, the stress am- 
plitude a is increased exponentially in time from 0.1 Pa 
to 1000 Pa. These data are thus the same as those pre- 
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FIG. 7: Viscoelastic moduli G' (•) and G" (o) as a function 
of time for a 4.7 wt % N5HU-toluene gel. The frequency of 
the stress oscillations is 1 Hz. For t < 480 s, the amplitude of 
the applied shear stress is increased from 0.1 Pa to 1000 Pa. 
For t > 480 s, a periodic stress of amplitude 1 Pa and fre- 
quency 1 Hz is applied. The gel recovers instantaneously its 
mechanical properties. 



sented in Fig. \5\ where they were shown as a function of 
the stress amplitude. The gel breaks at t ~ 300 s for 
cr ~ 60 Pa and then behaves like a viscous fluid (G" ~ 
for 300 < t < 480 s). 

For t > 480 s, a small periodic stress of amplitude 
1 Pa and frequency 1 Hz is applied. Instantaneously, 
the sample recovers its elastic behavior. This gel clearly 
exhibits a thixotropic behavior as defined in 1975 by the 
British Standards Institution Q: the apparent viscosity 
decrease under stress is followed by a gradual recovery 
when stress is removed l2ll. 



B. Global flow curve 

We now study the nonlinear behavior of our gel under 
constant applied shear rate. As previously, a 4.7 wt % 
solution of NII5U in toluene is prepared and heated at 
50° C in order to totally dissolve the solid powdered sam- 
ple. The hot solution is introduced between the two coax- 
ial cylinders of a Couette geometry made of PVC (gap 
e — 0.5 mm, inner radius Ri = 24.5 mm and height 
H = 30 mm) and then cooled down to 5°C with a cool- 
ing rate of about 20°C/min using an icy water circulation 
around the outer cylinder. The reason for using a Cou- 
ette geometry in this section is that we wish to further use 
ultrasonic vclocimctry (which is currently available only 
in Couette geometry) in order to construct a local flow 



FIG. 8: Global flow curve: shear stress o" as a function of the 
applied shear rate 7 for a 4.7 wt % N5HU-toluene gel. Two 
different regimes are evidenced below and above 7 ~ 70 s^^. 

curve by combining global rheological data and local ve- 
locity measurements. The same behaviour was observed 
in cone-and-plate geometry even though the global flow 
curve may vary a little bit. 

Using the shear rate imposed mode of the ARIOOON 
rheometer, a constant shear rate 7 is applied to the sam- 
ple for 10 min and the shear stress a{t) is recorded as a 
function of time. This procedure is applied for 7 ranging 
from 0.05 to 1400 s~^. Figure [5] presents the global flow 
curve (T versus 7 (also referred to as the "engineering" 
flow curve in the literature) where the measured shear 
stress values are obtained by averaging a{t) over the last 
300 s of each shear rate step. Waiting for 5 min be- 
fore starting to average the rheological data ensures that 
steady state is reached. Errors bars correspond to the 
standard deviation of a{t) and indicate noisy or oscilla- 
tory signals. 

The flow curve of Fig. [5] may easily be divided into two 
different regimes. First, for 0.05 < 7 < 70 s^^ (regime 
1), the measured shear stress decreases as a function of 
the applied shear rate. In the case of pure shear flow, a 
negative slope in the flow curve cannot be observed ex- 
perimentally since it corresponds to material instability. 
Thus the global flow curve of Fig. [S] strongly suggests the 
presence of inhomogcneous flows in regime 1. Second, 
for i ^ '^0 (regime 2), the measured shear stress 
increases as a function of the applied shear rate, consis- 
tently with a low-viscosity liquid behavior. Note that, if 
we repeat the experiments after shearing the system at a 
shear rate of 800 for 2 min and stopping the shear for 
5 min, we get the same rheological behaviour. This con- 
firms the thixotropic nature of our gel. To better analyze 
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FIG. 9: Experimental setup for ultrasonic velocimetry under 
shear in Couette geometry. The notations are defined in the 
text. 

this striking flow curve, we have performed local veloc- 
ity experiments simultaneously to nonlinear rheological 
measurements. 



C. Experimental setup for ultrasonic velocimetry 

Our rhco-ultrasonic setup allows the recording of ve- 
locity profiles across the gap of a Couette cell simulta- 
neously to rheological data with a spatial resolution of 
about 40 ^m. The time needed to record one velocity 
profile ranges between 0.02 and 20 s, depending upon 
the value of the applied shear rate. This setup has been 
described at length in Ref. Q and only its main charac- 
teristics are recalled here. It is built around the stress- 
controlled rheometer (TA Instruments ARIOOON) used so 
far and equipped with the Couette cell described above. 
The whole cell is surrounded by water whose tempera- 
ture is kept constant to within ±0.1°C. A piezo-polymer 
transducer of central frequency 36 MHz generating fo- 
cused ultrasonic pulses is immersed in water in front of 
the stator, as sketched in Fig. [51 The transducer is con- 
trolled by a pulser-rcccivcr unit (Panametrics 5900PR) 
which generates 220 V pulses with a rise time of about 
1 ns. The pulse repetition frequency Jprf is tunable 
from to 20 kHz. Backscattered (BS) signals are sam- 
pled at 500 MHz, stored on a PCI digitizer with 8 Mb 
on-board memory (Acqiris DP235), and later transferred 
to the host computer for post-processing. 

Under the assumption of single scattering, the signal 
received at time tus can be interpreted as interferences 
coming from scatterers located at position y = Coi„s/2, 
where Cq is the sound speed in the fluid and y is the dis- 
tance from the transducer along the ultrasonic beam (see 
Fig. [5]). When the sample is submitted to a shear flow, BS 
signals move along with the scatterers. Two successive 
pulses separated by a time interval r = 1/ /prf lead to 
two similar backscattered signals that are shifted in time. 
Figure [10] shows forty successive BS signals recorded in 
a sheared Newtonian latex suspension with a pulse fre- 




FIG. 10: BS signals corresponding to 40 pulses sent every 
T — 6.7 ms in a Newtonian latex suspension sheared at 7 = 
10 s^^. The normalized pressure p/po is coded in gray levels. 
The horizontal axis is the (fast) ultrasonic time tus and the 
vertical axis is the (slow) pulse time t. The gap lies between 
tua — 15.1 /is and tus — 16.0 ^s. 



qucncy repetition Jprf = 150 Hz. The amplitude of the 
BS signals is coded in gray levels and the vertical axis 
t corresponds to the time at which each pulse is sent. 
The ultrasonic pulses leave the stator and enter the gap 
at tus — 15.1 /is and the rotor position corresponds to 
tus — 16.0 ^s. The slopes of the traces left by the echoes 
in this two-dimensional {tus , t) diagram are inversely pro- 
portional to the local velocities. Thus, the signature of 
shear is rather clear: velocities increase from the stator, 
where the scatterers remain almost fixed, to the rotor. 
Moreover, the presence of moving echoes for t^s > 16 ^s 
is not surprising: these echoes simply correspond to scat- 
tering of the wave reflected on the rotor. One can sec that 
the whole range tus = 15.1-16.0 fis gets covered by some 
significant speckle signal in about 20 ms in that case. 

From these BS signals, one can relate the time shift 
St{tus) between two echoes corresponding to two succes- 
sive pulses and received at time tus to the displacement 
of the scatterers Sy{y) = Cq St{tus)/'2 at position y. The 
velocity profile u„(y) (i.e., the y component of the veloc- 
ity vector v) is then simply given by Vy(y) = Sy{y)/T. A 
cross-correlation algorithm is used to estimate the time 
shift St as a function of tus and the sound speed Cq is 
measured independently (see Ref. for more details). 
In order to get a non-zero projection of the velocity along 
the acoustic axis, ultrasonic pulses enter the gap with an 
angle 9 ~ 15° as shown in Fig. O 

Using the notations of Fig. [HI the projection of the 
velocity vector v along the y-axis in Couette geometry 
is the sum of the projections of the radial component 
and of the tangential component Vg of v: 

Vy = VgSinO + v,.cos6 . (1) 
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Assuming the flow field to be purely tangential, one gets: 



Ve{x) 



v{x) , 



(2) 



where x is the radial position given by a; e— y cos 9. The 
determination of v{x) thus requires a precise determina- 
tion of the angle 9, which is achieved through a careful 
calibration procedure using a Newtonian fluid 0] . If the 
radial component v,. is non-zero, the term cos 9 comes 
into play in Eq. ([T]), so that the definition of v{x) given 
in Eq. ^ leads to 



v{x) = Ve{x) + 



v^{x) 
tan 9 



(3) 



Therefore, the interpretation of the experimental velocity 
profiles v(x) may become problematic whenever the flow 
is not purely tangential. It should also be mentioned that 
our setup is insensitive to any contribution of the vertical 
component to the velocity vector {i.e., in the vorticity 
direction z), since its projection along the acoustic axis 
is always zero. 

In principle, a velocity profile could be obtained by 
cross-correlating only two succesive BS signals. How- 
ever, in practice the speckle amplitude is never uniform. 
Locally, destructive interferences or the absence of scat- 
terers may lead to signal levels too small to be well ana- 
lyzed. To recover a full velocity profile, some averaging is 
then performed over several successive cross-correlations. 
Typically, on a stationary flow, 100 cross correlations 
are required to get a standard deviation of the veloc- 
ity measurements of less than 6% 0]. In the following, 
we will present both time-averaged data and "instanta- 
neous" measurements i.e. data obtained using only two 
successive cross-correlation. 

Finally, thanks to their specific microstructure made of 
rigid aggregates whose size is comparable to the acoustic 
wavelength used for velocimetry, our NII5U gels naturally 
scatter ultrasonic waves in the single scattering regime. 
Thus it was not necessary to seed the gel with scatter- 
ers to perform velocity measurements. The next section 
displays and analyzes the time-averaged data recorded in 
the NII5U organogel. 



D. Time-averaged velocity profiles 

In this section, we focus on the time-averaged velocity 
profiles recorded using ultrasonic velocimetry once rheo- 
logical variables have reached steady state for each shear 
rate step along the flow curve of Fig. [S] To this aim the 
local velocity measurements v{x) were averaged over 6 
to 10 series of 1000 pulses (except for shear rates smaller 
than 0.3 s^^ where 2 to 5 series of 500 pulses were used 
in order to fit within the last 300 s of each shear rate 
step). 15 measurement points are obtained over the gap 
width e = 0.5 mm correponding to a spatial resolution of 
33 /im. The standard deviation of these measurements 
then reflects the occurence of temporal velocity fluctua- 
tions. 




FIG. 11: Time- averaged velocity profile v{x) as a function of 
the reduced coordinate x/e for 7 = 0.05 s~^. The horizon- 
tal bar on the left axis indicates the rotor velocity. No wall 
slip can be clearly evidenced and the measurements are very 
noisy. The profiles have been averaged on the last 100 s of 
the shear rate step at 7 = 0.05 s^^ once rheological variables 
have reached steady state. 



1. Low shear regime (regime 1) 

Figures [TT] and [T^] show the resulting time-averaged 
velocity profiles for shear rates ranging within the de- 
creasing part of the flow curve (regime 1). By definition 
of x, X = (cc = e resp.) corresponds to the rotor (sta- 
tor resp.). Three main features show up. We first note 
that for all applied shear rates, the velocity of the sample 
vanishes near the stator and reaches the velocity of the 
moving wall near the rotor. We conclude that no wall 
slip occurs for this gel. 

Second, the standard deviations of the velocity mea- 
surements are huge: they may reach 100% of the mea- 
sured value at low applied shear rate (see errors bars on 
Figs . [TT] and [T2)) . In some cases "instantaneous" negative 
values of the velocity are even recorded (see Fig. [TT|) . The 
standard deviation decreases as a function of the applied 
shear rate. Such large error bars clearly reflect the oc- 
curence of strong temporal velocity fluctuations (since 
error bars of less than 3% would be expected for a sta- 
tionary flow with the same averaging [9|). We will come 
back to this point in details in Sect. IIIIFI below. 

Third, in the bulk, three different zones (or "shear 
bands") can be clearly distinguished, (i) Close to the 
stator, the gel does not flow and behaves like a solid (see 
x/e G [0.55,1] in Fig. [TT] and x/e £ [0.75,1] in Fig.[T2l). 
(ii) In the middle of the gap (see x/e G [0.4,0.55] in 
Fig. [m x/e G [0.3,0.75] in Fig. [Hlja), x/e G [0.2,0.75] 
in Fig. [T^b). x/e G [0.1,0.75] in Vig.^c), and x/e G 



8 




FIG. 12: Time- averaged velocity profiles v{x) in regime 1 as 
a function of the reduced coordinate x/e. (a) 7 = 0.1 (o), 0.2 
(□), and 0.3 s"^ (>). (b) 7 = 0.4 (o), 0.5 (□), and 1 s"^ (>). 
(c) 7 = 0.2 (o), 5 (□), and 10 s~^ (>). (d) 7 = 30 (o) and 
70 (□). The horizontal bars on the left axis indicate the 
respective rotor velocities. No wall slip is evidenced. 



[0.05,0.75] in Fig. [121 (1)). the gel is flowing and bears a 
non-zero local shear rate. The value of this local shear 
rate increases as a function of the engineering applied 
shear rate. In the following, this zone will be called 
the "viscous zone." (iii) Finally, close to the rotor (see 
x/e £ [0,0.4] in Fig. [TTl x/e £ [0,0.3] in Fig. ^a), 
x/e e [0,0.2] in Fig. [IHb), x/e £ [0,0.1] in Fig. [HJc), 
and x/e £ [0,0.05] in Fig. [T2y d)). the gel remains un- 
sheared. The size of this solid-like region located near 
the rotor decreases as the engineering applied shear rate 
is increased. At the end of regime 1 (see 7 = 70 s^^ in 
Fig. HHd)), the solid-like state near the rotor disappears 
and a flowing liquid coexists with a solid-like region near 
the stator. 

Figure [13] displays the thickness of the three shear 
bands and the local shear rates in the viscous zone as a 
function of the engineering applied shear rate in regime 1. 
The local shear rate is extracted from the time-averaged 
velocity profiles using: 



-fix) ^ -r — ' 
ar 



(4) 



where r = + x. These data clearly reveal that the 
growth of the viscous zone occurs in the decreasing part 
of the engineering flow curve (i.e., for 7 < 70 s~^). This 
is not the case in clay gels where the growth of this region 
was shown to take place at a given shear stress @. We 
also note that the flow behavior of this viscous zone is 
rather complex. Indeed for a given engineering applied 
shear rate, the local shear rate may vary by a factor of ten 




FIG. 13: Analysis of time-averaged velocity profiles in 
regime 1. Top: thickness of the shear bands divided by the 
gap width e as a function of the engineering applied shear rate 
7. □ correspond to the solid-like state located near the stator, 
o to the middle zone, i.e., to the viscous zone, and ■ to the 
solid-like state located near the rotor. Bottom: local shear 
rates borne by the sample in the viscous zone as a function 
of the engineering applied shear rate 7. 



within the viscous zone. More precisely the local shear 
rate increases as a function of the reduced coordinate 
x/e: highest shear rates are evidenced near the solid-like 
state close to the stator. This result will be discussed in 
more details in the next section. 



2. High shear regime (regime 2) 

Time-averaged velocity profiles measured in regime 2 
are shown in Fig. [T4| In this regime, a flowing zone al- 
ways coexists with a solid-like region near the stator. In- 
creasing the applied shear rate increases the shear rate 
borne by the flowing zone, but almost does not change 
the size of the solid-like region near the stator. Error 
bars are stifl noticeable in regime 2 but their magnitudes 
are much less than those recorded at low shear rate. 



E. Local flow curve 

To go deeper into the analysis of the flow behavior, one 
may extract a local flow curve from the combined veloc- 
ity and rhcological measurements 0, d, H^. Indeed, the 
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FIG. 14: Time-averaged velocity profiles v{x) in regime 2 as 
a function of the reduced coordinate x/e. 7 = 100 (o), 200 
(>), 400 (<), and 800 s"^ (□). The horizontal bars on the left 
axis indicate the respective rotor velocities. No wall slip is 
evidenced. A solid-like region is detected near the stator even 
for the highest applied shear rates. 



FIG. 15: Local flow curve o-{x) versus 7(2;) inferred from 
rheological and local velocity measurements. Three different 
branches are evidenced: a solid-like state (□), a viscous state 
(o) in regime 1, and a low viscosity shear thinning fluid {+) 
in regime 2. The solid line corresponds to the shear thinning 
behavior a = A7" with A = 0.04 SI and n = 0.6. 



torque r imposed on the moving cylinder by the rheome- 
tcr yields the stress distribution cr(x) across the Couettc 
cell, 



and the velocity profile leads to the local shear rate 7(0;) 
according to eq. jl]). 

The resulting <7{x) versus 7(2;) data arc plotted in 
Fig. [T5l The first striking observation is that the local 
flow curve is multivalued: for shear stresses ranging be- 
tween 0.7 and 4 Pa, three different shear rates correspond 
to the same shear stress. The gel may thus adopt three 
different states for a given shear stress: (i) a solid-like 
state (□), (ii) a viscous state (o), and (iii) a low viscosity 
shear thinning state State (ii) is observed in regime 
1 and corresponds to the viscous zone in the velocity pro- 
files of Fig. [T2I The shear thinning local behavior (iii) is 
observed in regime 2 (see solid line in Fig. [T5|l and cor- 
responds to the flowing region in Fig. 1141 

Second, we note that flowing samples are observed for 
shear stresses below the upper limit of the linear regime 
CT ~ 3 Pa. This simply means that during the experi- 
ments under oscillating shear stress, only the solid-like 
state was probed. Moreover, states (ii) and (iii) are ob- 
tained after a fracturing process which probably requires 
a certain amount of strain and thus may not have enough 
time to develop under stress oscillations. 

Finally the most important result is the presence of 



a decreasing branch in the local flow curve (o and o in 
Fig. fTS]). In regime 1, the shear stress decreases from 7 
to 0.7 Pa while the applied shear rate is increased from 
0.05 to 70 s~^. Usually the apparent shear viscosity is 
defined as the ratio between a{x) and 7(0;). Decreasing 
branches in the flow curve are mechanicaly unstable. It 
is thus striking to get sound measurements in this zone of 
the flow curve. Let us recall that the velocity profiles dis- 
play very large errors bars in this range of applied shear 
rates, reflecting important temporal fluctuations. In or- 
der to get further insight on this decreasing branch, it is 
thus necessary to focus on the "instantaneous" velocity 
profiles and to study in details the temporal behaviour of 
the local velocities. The next section is devoted to this 
point. 



F. Analysis of backscattered ultrasonic signals and 
"instantaneous" velocity signals 

In order to better understand the decreasing part of 
the local flow curve, we go back to the backscattered 
ultrasonic signals recorded in regime 1 . Figure ll6r top) 
displays backscattered ultrasonic signals recorded in the 
4.7 wt % NH5U-toluene gel during the experiment of 
Fig. [8] at the very beginning of the shear rate step at 
7 ~ 0.05 s~-^. More precisely, shear is applied at i ~ 140 s 
and 100 successive signals with fpRF — 2 Hz are shown. 
Once shear is started, a flowing zone is nucleated in the 




FIG. 16: Top: BS pressure signals recorded in a 4.7 wt % 
NH5U-toluene gel under an applied shear rate 7 — 0.05 s~^ 
with T = 0.5 s. Arrows highlight discontinuities in the BS 
signals versus the pulse time t, which we interpret as the 
signature of fracture-like events. Bottom: "instantaneous" 
velocity extracted from the above data as a function of time 
at a:: = 0.33 mm (corresponding to tus = 15.70 fis). Arrows 
indicate the fracture events. 

middle of the gap after a few seconds and grows in time 
until steady state is achieved. The flow is thus clearly 
inhomogeneous . 

A second very striking feature of the BS signals at 
7 = 0.05 is that the traces in Fig. [TBTtop) present 
sudden discontinuities from one pulse to another, i.e., as 
a function of the pulse time t (sec arrows in Fig. I16p . 
For the sake of clarity, wc have also plotted the tempo- 
ral evolution of the "instantaneous" velocity at a given 
point in the gap {x ~ 0.29 mm, which corresponds to 
tus = 15.7 fj,s) in Fig. [rBr bottoml . Such "instanta- 
neous" measurements were obtained by cross-correlating 
two successive signals without further averaging cross- 
correlations. Although this procedure leads to rather 
noisy velocity signals, very fast variations of the "instan- 
taneous" velocity (which in this case becomes negative) 
are revealed, corresponding to the sudden discontinutics 
pointed out by arrows in Fig. fTHT top). Moreover the fact 
that these discontinuities do not span over the full range 
of ultrasonic times t.^s is an indication that "fast" events 
occur that are localized in both space and time. The time 
scale for such events is smaller than the time interval r 
between two pulses. Even if more quantitative experi- 
ments would be needed (in particular on much shorter 
time scales) in order to draw definite conclusions, wc in- 
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FIG. 17: Top: BS pressure signals recorded in a 4.7 wt % 
NH5U-toluene gel under an applied shear rate 7 = 0.1 s~^ 
with r = 0.25 ms. Arrows highlight discontinuities in the 
BS signals versus the pulse time t, which we interpret as 
the signature of fracture-like events. A typical zone of three- 
dimensional motion is surrounded by a white frame. Arrows 
indicate fracture-like events. Bottom: 'instantaneous" veloc- 
ity extracted from the above data as a function of time at 
X = 0.2 mm (corresponding to tus ~ 15.46 fis). Arrows indi- 
cate the fracture events and the dotted frame points to nega- 
tive velocity data on "slow" time scales, which are indicative 
of three-dimensional motion. 



terpret these events as fractures in the bulk of the fluid. 
Such a behavior persists when steady state is achieved. 
Figure [T7] shows that discontinuities in BS signals are also 
found for 7 = 0.1 s~^. 

The third important point to notice is that the local 
slopes of the BS signals present temporal variations on 
"slow" time scales (typically a few time intervals r be- 
tween two pulses). In some cases, the slopes do not 
only vary but also change sign even when steady state 
is reached. Such events are highlighted for 7 = 0.1 s~^ 
by white frames in Fig. [T71 They are also found at other 
shear rates in regime 1 (see t ~ 150-160 s in Fig. [TCTtop') 
and Fig. [TO]) . This apparent "reverse" motion of the ul- 
trasonic echoes can be attributed either to some rota- 
tional motion of the N5HU aggregates or to some signif- 
icant radial component of the velocity field. Although 
we cannot discriminate between these two mechanisms 
through the present experiments alone, our observations 
constitute a clear evidence for three-dimensional features 
of the displacement field with characteristic times of a 
few seconds. As explained above in Sect. IIII C[ the 
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15.2 15.4 15.6 15.8 16 



FIG. 18: BS pressure signals recorded in a 4.7 wt % NH5U- 
toluene gel under an applied shear rate 7 = 1 s^^ with r — 
50 ms (top) and 7 = 70 s~^ with r = 1.4 ms (bottom). 
Each subfigure displays a series of 20 successive signals. Some 
typical zones of three-dimensional motion are surrounded by 
a white frame. 



consequence of the three-dimensional nature of the flow 
is the possibility to measure apparent negative veloci- 
ties on rather long time scales. This is illustrated in 
Fig. llZr bottom) where the velocity remains negative for 
t ~ 92-94 s. 

Three-dimensional motions and fractures persist 
throughout regime 1 but these events become more rare 
and weaker as the shear rate is increased. Figure [TH] 
shows more examples of three-dimensional flow features 
recorded in regime 1 at 7 = 1 and 70 s^^. However no 
fracture-like events can be detected on these BS signals. 
As illustrated by Fig. [121 three-dimensional behaviors as 
well as fractures completely disappear in regime 2. 

This qualitative analysis allows us to better understand 
the huge errors bars reported in the time-averaged ve- 
locity profiles. As suspected, these error bars originate 
from very large temporal fluctuations in the flow field. 
We have shown that these fluctuations arise either from 
"fast" fracture-like events (on time scales smaller than r) 
or from "slow" three-dimensional motions (on time scales 
of a few r). Such features are not found at higher shear 
rates. We believe that such fluctuations cancel out when 
time-averaged data arc considered so that the velocity 
profiles of Figs. [TT] and zone It may still be interpreted in 
terms of tangential velocity. 

From these observations we may propose an explana- 
tion for the negative slope of the local flow curve reported 
in Sect. IIII El Indeed, fractures and three-dimensional 
flow features lead to enhanced friction losses in the sam- 
ple. Since their number and amplitude were shown to 
decrease as the shear rate is increased in regime 1, the 
loss of energy becomes smaller, which could lead to a de- 
crease of the shear stress. Note, however that as we are 
unable to submit the gel to pure shear flow, it is thus 
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FIG. 19: BS pressure signals recorded in a 4.7 wt % NH5U- 
toluene gel under an applied shear rate 7 — 200 s~^ with 
r — 0.25 ms. Each subfigure displays a series of 20 successive 
signals. No three-dimensional motions or fracture-like events 
are evidenced. 

not possible to identify a{x)/'j{x) with the shear viscos- 
ity for 7 < 70 s~^. We stress again the fact that the flow 
behavior along the decreasing branch of the flow curve is 
very complex since widely different local shear rates are 
evidenced in the viscous zone for the same value of the 
shear stress. 

Finally note that the exact location of the fractures or 
of the shear bands may depend upon the interactions of 
the gel with the walls. Experiments dealing with another 
orgaiiogel have shown that the fracture may locate at the 
wall [231. In this case, the gel was not cooled directly in 
the Couettc cell but before being loaded into the cell 
which clearly limits the links between the gel and the 
cell. 



IV. SUMMARY, DISCUSSION AND 
CONCLUSION 

In this work, we have studied the nonlinear rheolog- 
ical properties of an organogel. Under small applied 
shear rates (regime 1), a flowing region is nucleated in 
the middle of the gap of a Couctte cell. Increasing the 
applied shear rate increases the size of the flowing re- 
gion. Moreover backscattered ultrasonic signals clearly 
evidence three-dimensional displacements and fracture- 
like events in this "viscous zone." Time-averaged velocity 
measurements show that the local flow behavior in the 
viscous zone is characterized by a decreasing local flow 
curve and a very large range of accessible shear rates for a 
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given shear stress. At higher appUed shear rates (regime 
2), an increasing flow curve is recovered but a sohd-hke 
region persists close to the stator. When the mechani- 
cal stress is released, the fractures heal rapidly and the 
sample recovers its initial elastic properties. 

We believe that the present work sheds new light on 
the jamming concept and wish to further discuss our re- 
sults in this context. Indeed the jamming phase diagram 
suggested by Liu and Nagel [1^ provides a unified pic- 
ture of the ergodic-to-nonergodic transition in systems as 
diverse as cla y su spensions j2^ , emulsions [1^ , or micro- 
gel solutions [26[ . This diagram considers the possibility 
of unjamming the system by increasing its temperature, 
releasing its density, or increasing the amplitude of some 
applied shear stress. Although the jamming transition 
presents some universal features (universal slow dynam- 
ics and ageing behavior have been observed in a wide 
variety of jammed systems such as colloidal suspensions, 
concentrated emulsions, and lamellar gels), we argue that 
the unjamming mechanisms are not universal. In a pre- 
vious work, we have pointed out that the flow of concen- 
trated emulsions in the vicinity of the yield stress may be 
homogeneous or heteregeneous depending on the inter- 
actions between the droplets: shear banded flow occurs 
only for adhesive emulsions whereas the flow of nonadhe- 
sive emulsions remains homogeneous through the yielding 
transition 0. 

With the present work, we bring new elements that 
corroborate this point of view. First we have found that 
the viscous (or "unjammed") zone in regime 1 is nucle- 
ated in the middle of the gap of a Couette cell. Such a 
location for the flowing material had not been reported 
before: previous works have shown that, in the vicinity 
of the yield stress, flow is usually initiated in the region 
of higher shear stress i.e., close to the inner cylinder of 
a Couette cell Q. Moreover for shear rates in regime 1, 
the flow in the viscous zone is quite unusual with three- 
dimensional features and fracture-like events. Rotational 
motions of the aggregates and fractures induce a loss of 
energy. These events decrease in number and in ampli- 



tude as the applied shear rate is increased. Thus a lower 
shear stress is required to maintain the flow at higher 
shear rates. This scenario could provide an explanation 
for the decreasing branch observed in regime 1 on the lo- 
cal flow curve. We also point out that the jammed phase 
and the unjammed one may coexist over a very large 
range of shear stresses. In this case, defining and measur- 
ing a "yield stress" becomes very difficult. For shear rates 
in regime 2, fractures and three-dimensional motions dis- 
appear. Although the global flow curve points to a low 
viscosity liquid behavior, local measurements show that 
the system remains jammed close to the stator over a sig- 
nificant thickness and for the highest investigated shear 
rates. In the framework of the jamming phase diagram, 
this last feature remains rather puzzling. 

This study clearly shows that local measurements are 
required to understand the flow of complex fluids such 
as organogels. It also demonstrates that complex three- 
dimensional flows may occur in viscous systems, even at 
very low shear rates. In particular, although restricted 
to a specific organogel with large-scale heterogeneities, 
the present work reveals the existence of both "fast" and 
"slow" spatially localized events involved in the unjam- 
ming transition. Instead of claiming for any generality, 
we rather wish to emphasize that yielding mechanisms 
may well be widely different from one material to an- 
other. In our opinion, investigating more deeply the na- 
ture of the interplay between the complex local flow be- 
havior observed along the decreasing branch of the flow 
curve and the three-dimensional motions and fracture- 
like events revealed by the ultrasonic BS signals con- 
stitutes the most promising research direction on such 
thixotropic systems driven far from equilibrium. 
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